rbacon plot tests

require(rbacon)
#require(devtools)
#load_all("~/Dropbox/devsoftware/rbacon")
fry <- tempdir()
Bacon(run=FALSE, ask=F, coredir=fry, ssize=500, BCAD=F)

agedepth()

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm

## 100% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92

agedepth(rotate.axes=TRUE)

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm

## 100% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92

agedepth(BCAD=TRUE)

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm

## 100% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92

agedepth(rotate.axes=TRUE, BCAD=TRUE)

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm

## 100% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92

agedepth(rev.d=TRUE)

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm

## 100% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92

agedepth(rev.age=TRUE)

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm

## 100% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92

agedepth(rev.d=TRUE, rev.age=TRUE)

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm

## 100% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92

agedepth(rotate.axes=TRUE, rev.d=TRUE, rev.age=TRUE)

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm

## 100% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92

agedepth(BCAD=T)

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm

## 100% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92

agedepth(BCAD=T, rev.age=T)

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm

## 100% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92

agedepth(BCAD=F)

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm

## 100% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92

proxy.ghost()

## Calculating histograms

## 

layout(1)
accrate.depth.ghost()

accrate.depth.ghost(rotate.axes=T)

accrate.depth.ghost(rev.d=T)

accrate.depth.ghost(rev.acc=T)

accrate.depth.ghost(rev.d=T, rev.acc=T)

accrate.depth.ghost(rev.d=T, rev.acc=T, rotate.axes=T)

accrate.age.ghost()

## 

accrate.age.ghost(rotate.axes=T)

## 

accrate.age.ghost(rev.age=T)

## 

accrate.age.ghost(rev.acc=T)

## 

accrate.age.ghost(BCAD=T)

## 

accrate.age.ghost(BCAD=T, rev.acc=T, rev.age=T)

## 

accrate.age.ghost(rev.age=T, rev.acc=T)

## 

accrate.age.ghost(rev.age=T, rev.acc=T, rotate.axes=T)

## 

myflux <- cbind(0:100, rnorm(101, 100, 5))
flux.age.ghost(flux=myflux)

## 

## 

flux.age.ghost(rotate.axes=T, flux=myflux)

## 

## 

flux.age.ghost(rev.flux=T, flux=myflux)

## 

## 

flux.age.ghost(rev.age=T, flux=myflux)

## 

## 

Bacon(slump=c(20,21), ssize=500, coredir=fry, ask=F)

##  Will run 2,200,000 iterations and store 500

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 198 yr, min. 98 yr at 74.5 cm, max. 298 yr at 27.5 cm

## 100% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 23.26, shape 1.97, memory mean 0.52, strength 9.25

Bacon("RLGH3", coredir=fry, accept.suggestions=T, ask=F)

##  Using a mix of cal BP and calibrated C-14 dates

## Adapting acc.mean to 50 yr/cm

##  Will run 9,187,500 iterations and store 4,000

## Calculating age ranges...

## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 578 yr, min. 57 yr at 0 cm, max. 1835 yr at 220 cm

## Warning! Only 76% of the dates overlap with the age-depth model (95% ranges)

## Posteriors: accrate mean 46.08, shape 0.81, memory mean 0.47, strength 12.25

accrate.age.ghost(BCAD=T)

##